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Abstract: The Soil Moisture and Ocean Salinity (SMOS) mission is an Earth Explorer 
Opportunity mission from the European Space Agency (ESA). Its goal is to produce 
global maps of soil moisture and ocean salinity using the Microwave Imaging Radiometer 
by Aperture Synthesis (MIRAS). The purpose of the Passive Advanced Unit Synthetic 
Aperture (PAU-SA) instrument is to study and test some potential improvements that could 
eventually be implemented in future missions using interferometric radiometers such as the 
Geoestacionary Atmosferic Sounder (GAS), the Precipitation and All-weather Temperature 
and Humidity (PATH) and the Geostationary Interferometric Microwave Sounder (GIMS). 
Both MIRAS and PAU-SA are Y-shaped arrays with uniformly distributed antennas, but the 
receiver topology and the processing unit are quite different. The purpose of this work is 
to identify the elements in the MIRAS 's design susceptible of improvement and apply them 
in the PAU-SA instrument demonstrator, to test them in view of these future interferometric 
radiometer missions. 

Keywords: microware; interferometric radiometer; calibration; soil moisture and ocean 
salinity; SMOS; Passive Advanced Unit Synthetic Aperture (PAU-SA) 



Sensors 2012, 12 



7739 



1. Introduction 

Soil Moisture (SM) and Sea Surface Salinity (SSS) may seem to be unconnected, but both variables 
are intrinsically linked to the Earth's water cycle and the climate system. Thanks to the technological 
advances, the interest of the scientific community in remotely measuring SSS and SM has increased in 
the last years, spending much effort in this direction. The three current and planned Earth observation 
missions focused in these variables are: the MIRAS instrument aboard the SMOS mission of the 
European Space Agency (ESA) launched on November 2, 2009 [1], the Aquarius radiometer from 
National Aeronautics and Space Administration (NASA) aboard the Satelite de Aplicaciones Cientificas 
(SAC-D) spacecraft of the Comision Nacional de Actividades Espaciales (CONAE) [2] launched on 
June 10, 2011, and the Soil Moisture Active and Passive (SMAP) mission of NASA [3] to be launched 
in 2014. All three missions carry an L-band microwave radiometer as primary instrument. A microwave 
radiometer is an instrument that measures the spontaneous electromagnetic radiation emitted by all 
bodies at a physical temperature different from 0 Kelvin. Microwave radiometers were first used in 
radio-astronomy in the 1930s [4]. In [5] the Passive Advanced Unit for ocean monitoring (PAU) project 
was proposed aiming at demonstrate that collocated measurements of sea brightness temperature and 
reflected Global Navigation Satellite System Reflectrometry (GNSS-R) signals can result in a significant 
improvement of the retrieved SSS. Its scientific goals are to perform ocean monitoring by passive remote 
sensing to improve the knowledge of the relationship of the GNSS-R signal with the sea state, and to 
improve the knowledge of the relationship between the L-band brightness temperature and the sea state. 
To accomplish these goals, the PAU sensor consists of three instruments that operate in a synergetic way: 

• an L-band radiometer to measure the brightness temperature, 

• a reflectometer to measure the sea state using reflected Global Positioning System (GPS) 
opportunity signals, sharing the Radio Frequency (RF) front-end, with the radiometer, and 

• an Infrared Radiometer (IR) to measure the physical surface temperature. 

The PAU project is also a test-bed of new technological demonstrators such as real aperture 
radiometers with digital beamforming and polarization synthesis [7], fully-digital synthetic aperture 
radiometers, etc. To perform this, a ground-based instrument demonstrator PAU-SA has been designed, 
implemented, and tested to validate these possible improvements. This work is organized as follows: in 
order to understand better these improvements, some basic concepts on interferometric radiometry are 
reviewed in Section 2. The interested reader is referred to [8,9] for a detailed review. Section 3 provides 
an overview of the PAU-SA instrument, and describe in detail the instrument itself. Section 4 describes 
briefly the processing steps of the PAU-SA data. The interested reader is referred to [10] for more 
details. Section 5 presents an comparison between the PAU-SA and the MIRAS instruments to illustrate 
which improvements and new techniques have been implemented and tested, and their rationale. Finally, 
Section 6 presents some measurements performed, and Section 7 summarizes the main conclusions and 
the lessons learned. 
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2. Basic Concepts on Interferometric Radiometry 

The fundamental operation of a synthetic aperture radiometer is the complex cross-correlation of the 
signals collected by each pair of channels (antenna + receiver), or baseline as shown in Figure 1. 

Figure 1. Sketch of the fundamental operation of an interferometer: each pair of receiving 
channels and a complex correlator form a baseline that measures a sample of the visibility 
function from [6]. 
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The variables in Figure 1 are: F^ q n (^ rj) the normalized antenna patterns of channels m and n, at p 
and q polarizations, (£, rj) = (sin 6 cos 0, sin 6 sin 0) are the direction cosines defined with respect to the 
X and Y axes, TA m ,n is the antenna temperature, H m , n (f) is the frequency response, G m?n is the power 
gain of the channel, £? m?n is the equivalent noise bandwidth, UP£ is the analytical signal collected by 
each channel, and () is the expectation operator implemented in practice as a time average. A synthetic 
aperture radiometer measures the cross-correlations between all pair of signals, collected by the antennas. 
The auto-correlation and the cross-correlation of the signals at each receiver output are defined as: 

\{K,n(t)\ 2 ) 4 k B ■ B m , n ■ G m , n ■ (TZ, n + T RECm , n ) (1) 

and 

\(P m (t)bl*(t)) ± k B • y/B m • B n • ^jG m • G n • V%(u mn , v mn ) (2) 

where k B is the Boltzmann constant, T REC m, n is the equivalent receiver's noise temperature, and V£? n is 
the visibility function defined in the spatial frequency (baseline) (u mn , v mn ) = (x n — x m , y n — y m )/Xo 
that depends on the difference of the antenna positions normalized to the wavelength, being A 0 = c/f 0 . 
According to Equation (2), the visibility function with units of Kelvin can be derived as a sample of the 
cross-correlation function: 

VZ(Umn,V mn ) = -— i— — — l(6^t)6;(t)) (3) 

Assuming that the collected signals are stationary random processes (thermal noise radiation) that 
fulfill the ergodicity property, narrow-band, spatially uncorrected, and that the distance between 
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antennas is much larger than the wavelength, Equation (3) can be computed in practice using the 
cross-correlation (time average) of the in-phase and quadrature components Q) p r ^ n = I%{ q n (t) +jQ p m q n {t))\ 
V™ n {u mn ,v mn ) = < IUI > +j < QUI >. 

From the Van-Cittert Zernicke theorem [11], the visibility function can be related to the brightness 
temperature distribution [12] as: 

2 n n Vj) — T rec 8qp 



vz(u mn ,v mn ) = -== // ;; " V ^v)Fr(^v) 

\J\L m \L n JJp+nZKi Wl — q — V 



fo 



where 



• Q m ,n are the equivalent solid angles of the antennas, 

• v) i s ^e T B of the scene at p — q polarization [13] (E p and E q being the electric fields at p 
and q polarizations), 

• T rec is the physical temperature of the receiver (the so-called Corbella's term) [12], 

• S pq is the Kronecker's delta function: 5 pq = 1 if p = q and S pq = 0 if p ^ q, 

Umni + VmnV | j s ^ S0 -called fringe-washing function. This term is related to the limited 

Jo J 

bandwidth and the frequency response of the filters in the two receivers forming the baseline, being 

r mn ± JZ°H m (f)Hn(f)*e> 2 **d f , and 

• 1 / y/l — £ 2 — rf is the obliquity factor. 

In addition, Equation (4) is normalized by the term k B • \ZB m B n • ^G m G n • y/Q m Q n to ensure that 
the visibility function has units of Kelvin. 

2.7. Ideal Situations 

A simplified version of Equation (4) is analyzed. Given the narrow bandwidth of PAU-SA (2.2 MHz) 
and the reduced dimensions of the array (being Ar & 2.4 m the maximum distance between two antennas 
of the array), the fringe- washing function becomes: 

7 I ~ 1 (5) 

JO J 

Assuming that all the antenna patterns are identical F m = F n9 and the constant terms are normalized to 
the unity, the relationship between the visibility function V mn , (Equation (7)) and the so-called modified 
brightness temperature (Equation (6)) is given by Equation (8): 

T(C,V) = ^M== 2 \Fn(tv)\ 2 (6) 



/ t F n (tv)\ e- j2v(Umni+Vmn,1 > d£dr) (7) 

In these conditions, the modified brightness temperature can be recovered by means of the inverse 
Fourier transform of the visibility samples. 

nt ) r l )=&- l [V mn ] (8) 
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3. PAU-SA Description 

3.1. PAU-SA Instrument Overview 

Figure 2 shows the PAU-SA instrument without the radome. It is composed by a Y-shaped array 
of 25 antennas for radiometric applications: 8 antennas uniformly spaced every d = 0.816 A per arm 
plus the one in the center. An additional antenna connected to a matched load, but without receiver or 
antenna dummy is included at the end of each arm to improve the antenna pattern similarity. Moreover, 
in the original design, the reflectrometer part (PAU-SA GNSS-R) was going to use the 4 central antennas 
plus 3 additional ones (7 antennas in total) to create a steerable array to point to the Global Navigation 
Satellite System (GNSS) signal specular reflection points. Nowadays, this part has been simplified by a 
single additional antenna plus a GPS receiver front-end and data logger for off-line processing. 



Each dual-polarization RF front-end is integrated behind a dual polarization patch antenna. The three 
main features of this receiver are: simultaneous dual polarization (V & H), frequency of operation (GPS 
LI band) the same for both instruments (radiometer, and reflectometer), a three-stage down-converter, 
and a Local Oscillator (LO) inside each down-converter generated from a 10 MHz reference master 
clock common to all receivers. In this way, since the LO is generated inside each down-converter, 
their noises are uncorrected from receiver to receiver and do not introduce a correlation offset. Each 
receiver translates the input signal from 1,575.42 MHz to an intermediate frequency of 4.309 MHz with a 
bandwidth of 2.2 MHz and an approximated gain of 1 10 dB. The differential Intermediate Frequency (IF) 
signal is sent to the Analog to Digital Converter (ADC) array unit through twisted pair (RJ45 grade 5) 
cables. These signals are digitalized at eight bits using IF sub- sampling techniques with a sample 
frequency of F s = 5.757 MHz. These data enters into a Virtex 4 Field Programmable Gate Array (FPGA) 
for digital In-phase and Quadrature (I/Q) down-conversion and digital filtering using 8 bits. Then, all 
possible cross-correlations between antenna pairs at each polarization are simultaneously performed at 



Figure 2. Picture of PAU-SA's array without radome from [6]. 
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1 bit only, using the sign bit, with maximum integration times multiple of 1 s. At the same time the 
8 bit digitalized signals are squared and accumulated to compute the signal power, instead of using 
power detector diodes, as in MIRAS. Using digital techniques it is possible to eliminate 3 possible error 
sources: quadrature errors in the I/Q analog demodulators, thermal drifts in the Power Measurement 
System (PMS) due to the diodes' thermal stability used in analog systems, and quasi-perfect matching 
of the frequency responses for all receiver's using a digital filter (the only differences are due to the 
RF analog filters, which are wider). Once the correlation matrix is calculated, the data is sent via a 
Universal Asynchronous Receiver/Transmitter (UART) interface from the FPGA to the internal Personal 
Computer (PC) for an off-line post-processing in the external PC Figure 3. Due to the high volume 
of data, it is not possible to process it inside the FPGA. The external PC performs three tasks: 
(1) to undertake the calibration process needed before each acquisition process; (2) to do the off-line 
post-processing for the image reconstruction; and (3) to record the data for off-line post-processing. 



Figure 3. Global view of the PAU-SA's architecture from [6]. 




The detailed description of the calibration approach is described in [10]. For calibration purposes, 
a noise distribution network distributes correlated noise into all receivers at two different noise levels 
(hot and warm) for phase calibration. Moreover, an uncorrelated load is used inside each receiver for 
correlation offset estimation. 

3.2. PAU-SA Instrument Description 

An overview of the PAU-SA instrument is first given in Section 3 to provide a global understanding 
of the system. Figure 4 shows the main block diagram where the main modules are interconnected, 
which will be described in detail in the next subsections. These modules have been separated in two 
main blocks: the PAU-SA instrument and the mobile unit, adopting a centralized control through the 
external computer. This basically makes two control tasks: to control the motion of the articulated 
arm, and to control the internal computer that at the same time controls the devices located inside the 
PAU-SA instrument. 
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Figure 4. PAU-SA's system block diagram indicating the interface connections between 
different modules from [6]. 
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3.2.1. Antenna Array 

The first element of the instrument are the antennas. These are square patch antennas [14] with 
dielectric air to resonate at 1,575.42 MHz (GPS LI frequency). The patch has been printed on 0.6 mm 
FR4 substrate offering simplicity and low production costs. The main parameters are: matching better 
than —22 dB in each polarization (V/H), low ohmic losses t]q = 0.98, and half-power beam width 
of 60°. The antenna array consists of 25 antennas distributed over a Y-shaped array, with 8 antennas per 
arm, and a central one for radiometry applications. Three extra dummy antennas (with no receiver 
connected) are placed at the end of each arm to improve the antenna voltage patterns similarity as 
analyzed in [15] and [16]. Moreover, three antennas are placed around the central one for future GNSS-R 
applications [17] (Figure 2). Nowadays, this part has been replaced by a ceramic patch GPS Left Hand 
Circular Polarization (LHCP) antenna and a GPS LI receiver [18] connected to the internal PC for 
acquisition of the raw GPS data. A maximum distance of X/^ (d = 0.577 A) between antennas is 
required to avoid aliasing in the image reconstruction. However, due to the physical dimensions of 
the receivers and antennas, the minimum feasible distance in PAU-SA is 15.5 cm (d = 0.816 A), thus 
some aliasing will affect the borders of the recovered image. 
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3.2.2. Receiver 



The receiver has been designed to combine two different sub-systems: an L-band radiometer and 
a GNSS-Reflectometer. Since the reflectometer part requires continuous data acquisition, the input 
cannot be chopped. Therefore, Total Power Radiometer (TPR) topology for each vertical and horizontal 
polarization has been selected (Figure 5). Each chain is individually amplified, down-converted, and 
amplified again to allow an 8-bit quantification for later digital processing (cross-correlation, power 
detection, and calibration). Since the noise introduced by the amplifiers is uncorrelated (not considering 
the cross-talk), it vanishes when the output signals are cross-correlated. 

Figure 5. PAU-SA's receiver block diagram uses two TPR topologies, one per polarization 
(V&H)from [6]. 
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To improve the Alias-Free Field Of View (AF-FOV), the receiver has been implemented in reduced 
dimensions (11 cm x 7 cm x 3 cm). Taking into account that the operating frequency is the LI 
GPS band, many already existing commercial off-the-shelf (COTS) components have been used for the 
implementation of this receiver, such as the Zarlink GP2015 GPS front-end as down-converter [19], Low 
Noise Amplifier (LNA), and Surface Acoustic Wave (SAW) filters. The receiver was carefully designed 
to preserve symmetry, minimize cross-talk, and interconnection routes, while maintaining equal delays 
for all paths, being divided in two parts: RF and IF stages. Due to space limitations the two stages have 
been implemented in two boards one over the other, using different substrates according to the frequency 
of each stage. 
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3.2.2.1. RF Stage 

In order to ensure a constant impedance of the transmission lines, the RF stage has been implemented 
with ROGERS 4003 0.8 mm thickness substrate. Its structure consists of two parts, the switching stage 
and then an amplification stage. The main objective of the RF stage is to achieve a gain of at least 
30 dB over the input signals, keeping the Noise Figure (NF) as low as possible. This is necessary 
to work in the linear region of the down-converter module located in the IF stage. The matching of 
each input/output ports have to be at least —22 dB. The required isolation between antenna signal and 
calibration (correlated and uncorrelated) signals has been experimentally determined to be better than 
— 80 dB. Finally, the cross-talk between adjacent TPR have to be better than —40 dB. 

The switching stage is controlled through CTR1 and CTR2 signals by the Control Unit (CU) 
implemented in the FPGA. It switches between three different states: the antenna measurements, 
calibration with external correlated signals (either noise or Pseudo-Random Noise (PRN) signal), or 
calibration with uncorrelated noise generated by an internal 50 Q matched load. Usually this stage is 
connected to antenna acquisition, necessary for GNSS-R applications. The antenna input is connected 
directly to the patch antenna through a 50 Q cable. Calibration with correlation signal is required to 
compensate different phases and amplitudes among channels. To implement this structure two switches 
model RSW-2-25P by MiniCircuits have been used, one for each radiometer, increasing the isolation 
between correlated noise and antenna signal. Once selected, the signal is amplified to adjust the antenna 
input power (~ — 1 10 dBm) to the down-converter linear behavior power margin. The total RF gain 
is ^33.3 dB, obtaining an input power level for down-converter of ^—76.7 dBm, within the linear 
behaviour, and the total NF of 2.2 dB. 

3.2.2.2. IF Stage 

The main goal of the IF stage is to shift the RF signal to IF using a superheterodyne receiver or 
down-converter. Moreover this stage pre-amplifies the IF signal and transmits it to the ADC module 
through an Ethernet cable. This stage has been implemented with FR4 substrate. The output of each RF 
chain is interconnected by means of a semi-flexible cable to the input of each down-converter located 
in this stage. It mainly consists of the translation from 1,575.42 MHz to 4.309 MHz with a bandwidth 
of 2.2 MHz amplifying the RF signal by approximately 52 dB. The down-converter requires a 10 MHz 
Transistor-Transistor Logic (TTL) clock signal. This signal is distributed by a coaxial cable and it is 
used to internally synthesize the different LOs by means of the Phase Lock Loop (PLL)s. Usually these 
devices provide the output at IF, digitized at 2 bits (sign and magnitude) enough to recover the Coarse 
Acquisition (C/A) code, used in GPS applications. Originally, it was intended to determine the signal's 
power using just 1 bit as in [20] using the Nemerix NM1 100, but this coarse sampling was found totally 
unacceptable for radiometric applications. Therefore, the Zarlink GP2015 GPS down-converter was used 
instead, since it has an analog output test that can be sampled. A commercial GPS down-converter can 
be used for radiometric applications if it has a linear behavior. However, all commercial GPS receivers 
have an Automatic Gain Control (AGC) circuitry that modulates the gain depending on the input signal 
to achieve a constant power level at its output. The GP2015 AGC control has been tuned forcing a 
differential voltage between its pins. In this case, for input levels from —85 dBm to —60 dBm it behaves 
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linearly, with a gain of 52 dB, and with a maximum linearity error of 0.25 dB (Note that the dynamic 
range of the radiometer is actually much smaller (less than 3 dB) for which the linearity 0.2%). 

At the receiver's output the signals are centered around f IF = 4.309 MHz, with a B = 2.2 MHz 
bandwidth. Each down-converter output is then amplified using a NE592D video amplifier, with fixed 
gain and differential output in order to cancel the common mode errors in the transmission to the ADCs 
array using Ethernet category 5 cables. Due to component tolerances it is necessary to perform manually 
a relative adjustment of each channel by means of the voltage divider, between the down-converter, and 
the video amplifier (Figure 4). The goal of the voltage divider is to adjust the gain to obtain a 110 mV 
rms signal at the receiver's output, to match the ADC input range. The gain of each stage should have a 
maximum error of ±1 dB to allow the correction in the digital process part. Due to the space limitations, 
the correlation function is implemented digitally outside the receiver, using a Virtex 4 LX60 FPGA [21]. 
The analog signals are sent to an external ADCs array to be digitalized. To do this, the output of each 
video amplifier is connected to a matching network to match its output impedance to the 150 Q of the 
twisted pair (RJ45 grade 5 cable). This connector is shared by four differential outputs at 4.309 MHz, 
two RF switch control inputs, and two power supplies. Due to the amount of pairs, two RJ45 connectors 
are needed. This method minimizes considerably the total number of cables connected to the receiver, 
and at the same time provides a high isolation against possible interferences. 

3.2.2.3. PAU-SA's Receiver Implementation 

The design of PAU-SA's receiver has been based on PAU-Real Aperture (PAU-RA) [7,22], the real 
aperture version of the PAU instruments family. However, due to the large number of receivers needed in 
the PAU-SA instrument and the limited number of I/O pins in the FPGA to implement the correlations, 
it was not possible to use the same pseudo-correlation topology, and a total power radiometer topology 
was selected instead of a pseudo-correlation one. As it can be seen in Figure 6(a) and 6(b), PAU-SA's 
receiver has only half of the components mounted in the IF and RF stages. 

Figure 6. (a) PAU-SA's receiver IF stage view with box, and (b) RF stage view without box, 
(c) assembled stages to minimize the required area, and (d) interconnection stages (from [6]). 
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Figure 6. Cont. 




Once the RF and IF boards are assembled separately, they have been assembled and interconnected, 
as shown in Figure 6(c) and 6(d). 

3.2.3. Digital Sub-Systems 

The sub-systems implemented in the FPGA are: 

• In-phase (I) and quadrature (Q) demodulation of the receivers' output digital signals coming from 
the array of 8 bits ADC), 

• Digital Low Pass Filter (LPF) at 8 bits, 

• Power estimation system of the 50 receivers' output signals (25 receivers x 2 polarizations @ 
8 bits). 

• Correlation unit of the three correlation matrices (V, H, VH) @ 1 bit, and 

• Communication protocol and control with a PC, and data collection. 

The FPGA device used in the design is the Xilinx Virtex-4 LX 60 distributed. Figure 7 shows the block 
diagram of the main sub-systems previously discussed and their implementation in the FPGA. Other 
peripherals connected to the FPGA, such as an external PC, and an array of ADCs are also shown for 
clarity. Before start explaining the design of the sub-systems implemented in a FPGA, it is necessary to 
analyze the analog-to-digital converter. Since there are 25 receivers, and each one has a dual polarization 
antenna with IF differential output, a total of 50 analog differential channels are connected to the array of 
ADCs. Moreover the receivers need two unipolar common control signals (CTRL1 and CTRL2) to select 
the different RF inputs. These two control signals are sent by the external PC to the receivers through the 
FPGA first, and then via the ADC array. The IF output signals are down-converted at a central frequency 
of 4.309 MHz with a bandwidth of 2.2 MHz, and a dynamic range of the output voltage of 1 V rms . 
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Figure 7. Overview of the sub- systems implemented in the FPGA (radiometer part) and 
peripherals in PAU-SA from [6]. 
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3.2.3.1. Sampling Frequency 

The fundamental requirement of a sequential electronic block is the frequency of operation. In the 
case of the ADC, the frequency of operation determines the sampling frequency. In the PAU radiometer 
the receiver has a 4.309 MHz intermediate frequency, and then it is digitized to translate the frequency 
to baseband with the Digital Down Converter (DDC). The Nyquist sampling theorem states that the 
sampling frequency should be at least twice the highest frequency components in the signal. This means 
that to sample a 4.309 MHz signal, a minimum 8.618 MHz sampling frequency is required. A higher 
sampling frequency implies a higher data rate, and the DDC function will be more complex to implement 
in the FPGA because the input signal of DDC should be multiplied by a cosine function using all the 
8 bits. Another option is the so called "band-pass sampling". With this technique a smaller sampling 
frequency can be used, thus reducing the data rate, while simplifying the implementation of the DDC. 
The "band-pass sampling" technique is described in the following paragraphs. If the input signal of 
DDC is multiplied by cos(^n), where Q = 2nF digita i is the carrier frequency, and F digita i is the digital 
frequency with a value of 0.25 (Figure 8(a)), the input signal results: 

cos(^ • n) = 1, 0, -1, 0, . . . with n = 0, 1, 2, 3, . . . (9) 



Figure 8. (a) Frequency domain of input signal, and (b) frequency domain of input signal 
digitalized. Pictures from [6]. 



2,2 MHz 




4.309 



(a) 



5 4 Frequency (MHz) 















{ -3dB 

-J 1 




\ I 


1 

0 0.25 


1 

0.5 


1 | 
0.75 1 


F digital 



V V 

(b) 



This means that now, the input signal has to be multiplied by 1,0,-1,0 thus simplifying the numerical 
computation. Now, it is possible to compute the sample frequency of the ADC, considering the sampling 
frequency of ADC (F s ) and the frequency of output signal from receiver (F IF ), the resultant frequency is: 



Ffinal — Fs — Fif 
Taking F FINAL = the F s results: 

F s — 4 • F final 4 • {Fs — Fif) F s 
Fif = 4.309 MHz F s = 5.745 MHz 



"IF 



(10) 

(11) 
(12) 



Thus, the output signal of the ADC is centered at the frequency F F inal = 1.436 MHz, as shown in 
Figure 8(b). 
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3.2.3.2. I/Q Demodulation Unit 

Before calculating the correlation matrices, it is necessary to down-convert the digital signals of the 
receivers to baseband and obtain their in-phase (I) and quadrature (Q) components. This section briefly 
presents the theoretical formulation and then explains the main blocks. The I/Q demodulation unit 
is composed of four blocks as shown in Figure 9: serial-to-parallel converter, two's complement, I/Q 
demodulation, and a selective LPF. The serial to parallel block is necessary to convert serial output of 
the ADCs forced by the limited number of pins in the FPGA to parallel for easier internal operations. 
Therefore the ADC device has been chosen with serial output and a sampling frequency of 5.745 MHz. 
Since the raw data that comes from the ADC is not suitable for data processing due its offset (128), 
this offset has to be subtracted. This offset is related to the fact that the ADC output gives a binary 
quantification of the analog signal between 0 and 255 (8 bits). Hence, it is necessary to subtract this offset 
to have symmetrical positive and negative values at the same time, convert the data into 2's complement 
to simplify the signal processing. The two main blocks are I/Q demodulation and the LPF, explained in 
detail in this section. 



Figure 9. Block I/Q demodulation unit from [6]. 
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Once the receiver's signal has been digitized the output has the following expression: 

S(n) = R{(i(n) + jq(n)) • e jQn } (13) 
where Q = 2nF digita i is the carrier frequency. The above given expression can be rewritten as: 

S{n) = i(n) cos(fin) — q(n) sin(fin) (14) 

Taking into a count the previous equation, it is possible to obtain the I(n) and Q(n) components 
multiplying S (n) by a cosine and a sine respectively, and low-pass filtering them: 



I(n) = S(n) cos(^n) = -i(n) H — (z(n) cos(2f2n) — q(n) sin(2f2n)) 
Q(n) = S(ri)(— sin(fin)) = -q(n) — -(i(n) sin(2fin) + q{n) cos(2fin)) 



(15) 
(16) 



If these expressions are particularized to the digital frequency (F digita i = 0.25), the product of S(n) 
by a cosine and a sine is equivalent to multiply S(n) by a periodical sequence with period four: 



,7T 



cos(2tt F d ig ital ri) = cos(-n) = 1, 0, -1, 0, . . . with n = 0, 1, 2, 3, . . . (17) 



,7T 



■ sm(27r Fdigitain) = -sin(-n) = 0, -1,0, 1, . . . with n = 0, 1, 2, 3, . . . (18) 
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Working with the ADCs and using band-pass sampling, it is then possible to work with a 
specific digital frequency that minimizes the required hardware resources. In this particular case 
it is not necessary to implement a multiplier stage, but only to take a sample out of two with the 
corresponding sign. These modules are implemented with a multiplexer block and an inverter function. In 
Equations (15) and (16) it is possible to observe a low frequency in-phase and quadrature terms plus a 
high frequency contribution that must be eliminated by a LPF. The filter implemented is an Infinite 
Impulse Response (IIR) LPF with a cut-off frequency of 0.25 in the digital domain with an attenuation 
larger than 20 dB at the frequency of interest, and 6 dB of gain that compensates the 0.5 factor in the 
phase and quadrature Equations (15) and (16), (Figure 10(a)). The disadvantage is that the phase does not 
have a perfectly linear behavior (Figure 10(b)), and the group delay changes with frequency. However, 
this is not relevant, since all filters have exactly the same behavior, and uniformity is guaranteed across 
them, unlike with analog filters. These filters have been implemented only with elementary functions: 
delay blocks, adders and shift registers using the minimum FPGA resources. 

Figure 10. (a) Gain response of the IIR filter (b) Phase response of the IIR filter. Pictures 
from [6]. 
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3.2.4. Power Measurement 

To recover the normalized visibility function as derived from the correlation matrix (Equation (2)), an 
estimation of the system temperature is necessary for both polarizations. This means 25 (receivers) x 2 
(polarizations) = 100 different power measures. To obtain this, each channel is considered as a total 
power radiometer and the signal's power is estimated as (Equation (19)) 

n 

Since our signal is thermal noise, the power measurement can be obtained either from the IF signal 
or from each of the I/Q components. For the sake of simplicity, the I and Q components have been used. 
To implement a power measurement system a multiplier and an adder are required. For this application 
the internally implemented FPGA functions are used. 
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3.2.5. Digital Correlation Unit (DCU) 

The computation of the digital correlations and the power measurements are the two pre-processing 
steps implemented in the FPGA. The Digital Correlation Unit DCU measures the similarity between both 
signals using 1 bit (sign bit). It basically consists of counting the number of samples (N c ) with the same 
sign between all possible baseline combinations. Moreover, it is also necessary to obtain the maximum 
number of samples (N cmax ), which is proportional to integration time, so as to normalize and obtain the 
correlation matrix (c) in the post-processing. Since the sampling frequency (F s ) is 5.745 MHz and that 
the maximum integration time (T int ) is 1 s, then the maximum number of counts (N cmax ) is given by 
Equation (26), requiring a 23 bits a 2 23 ~ 8.388 MS counter enough to count up to 5.745 MSamples. 
In this case, with 25 receivers and the I and Q demodulated components, the correlation matrix can be 
obtained (Figure 11). Recall that the correlations between the in-phase components (I-I) are the real part 
of the normalized visibility (// r ), and the correlations between in-phase and quadrature components (I-Q) 
correspond with the imaginary part of the normalized visibility (/x^). In the diagonal, the cross-correlation 
between the in-phase and quadrature components (I-Q) of the same signal are computed, but as expected 
due to the digital I/Q demodulation, they are zero. Basically, to implement the correlator block, two 
elementary components are required: a logic equality detector or XNOR gate, and a 23 bits counter, as 
shown in Figure 12. 



Figure 11. PAU-SA's correlation counts matrix AT Cm>n from [6]. 
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Figure 12. Elemental correlator block (from [6]). 
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To implement only one polarization matrix, 676 results of 23 bits are necessary. Due to limitations in 
the FPGA resources the internal Random Access Memory (RAM) memory has been used to implement 
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the counters (676 x 3 = 2,026 counters). Also, since the frequency of operation is higher than the 
sampling frequency (F C lk = 103.41 MHz 3> F s = 5.5745 MHz), it is possible to divide the sampling 
time in different slots, and reuse the hardware by a factor: 



T 

J. ry 



- clok 



18 



(20) 



This value means that it is possible to do 18 internal operations until the next sample arrive. In our 
case the FPGA used has 160 independent RAM blocks to implement 3 correlation matrices, this means 
53 blocks for each block. In this case, the required hardware reuse factor is: 

676 counters each matrix 



53 RAM blocks each matrix 
which is smaller than the maximum one (Equation (20)). 



12.7 



(21) 



3.2.6. Correlated Noise Unit 

Figure 13 shows the correlated noise diagram. It consist of three main blocks: two blocks for the 
signal generation (on the left), composed of thermal and pseudo-random noises, and a selection circuitry 
(on the right). Regarding the noise signals there are two different: a classical noise source, and a 
new technique based on using Pseudo-Random Noise PRN sequences, which can be selected using 
the selector block. 

Figure 13. Correlated noise unit block diagram from [6]. 
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The Classical noise source has been implemented with a NoiseCom noise source model NC346D with 
an Excess Noise Ratio (ENR) of 21.31 dB [23]. The ENR can be related with an equivalent temperature 
T n at the output of the noise source as: 



ENRdB^lOlog 



T 

-J- n 



1 [dB] 



(22) 
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where T 0 is the reference temperature (290 K). In this case, the resulting T n : 

1 J = 39.5 • 10 3 [K] 



EN Rj d 

T n = To • ( 10^- 



(23) 



3.2.6.2. Pseudo-Random Noise (PRN) 



Pseudo-Random Noise (PRN) sequences are signals with very long repetition periods that are used 
in a variety of applications, such as Code Division Multiple Access (CDMA) communications or 
positioning systems. They have a relatively flat spectrum over a bandwidth determined by the length of 
the sequence and the speed of the code or Symbol Rate (SR). Their spectra look like the noise spectrum, 
and the calibration of a microwave correlation radiometer (either interferometric or polarimetric) can 
benefit from these properties The SR parameter is used to determine the speed of the PRN code in order 
to control the bandwidth of the spectrum. The symbol rate is defined as the ratio of the bandwidth of 
the PRN signal (B PRN ) and the receiver's low-pass equivalent bandwidth (£?, Equation (24)). B PRN is 
related to the sequence duration (r PRN ) and the number of chips (a chip is like a bit, but it does not carry 
any information) N chips as shown in Equation (25). 



SR = 



B 



PRN 



B 



PRN 



B 

N 1 chips 
TpRN 



(24) 
(25) 



The equivalent noise temperature of the PRN signal (T PRN ) at the PRN generator module output is 
defined in terms of the PRN signal's amplitude (A): Pprn = 

A 2 /2 4 k B • T PRN • B PRN , where P P rn 

is the PRN signal power, and k B is the Boltzmann constant. 

The PRN is generated with a Linear Feedback Shift Register (LFSR) [24] as used for example in GPS 
applications. This part has been implemented in the FPGA model Spartan 3 of Xilinx using the module 
HLP-HS-FPGA [25] (Figure 14). 
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Figure 14. PRN generator module circuitry from [6]. 
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The system has been designed to have the possibility to select between 10 or 20 order primitive 
polynomials which maximum-length to generate two pseudo-random sequences of length 2 nblts — 1, 
respectively. The LFSR can produce an output at different speeds by selecting the SR parameter. 
Internally the reference clock of 1.023 MHz is multiplied by the SR parameter to change the speed 
of the sequence. Taking as a reference the polynomial of order 10 with a length sequence of 1,023 
chips, with SR = 1, the complete sequence is generated in 1,023 chips/1.023 MHz = 1 ms, with 
SR = 2, there are two complete sequences generated in the same time, and so on. With the 20th 
order polynomial with a length sequence of 1,048,575 chips, the sequence looks like more random, 
that is to say, with SR = 1, 1,048,575 chips/1.023 MHz = 1,025 ms is needed to repeat the sequence, 
and with SR = 2 half the time is needed, and so on. The output of the LFSR is baseband and it is 
internally modulated up to 80 MHz to be up-converted with the external mixer up to 1,575.42 MHz. 
Since the output of the digital signal is around 3.3 V, it has been necessary to put some attenuators in 
order not to saturate the receivers. In this case a 60 dB attenuator has been achieved cascading two 
attenuators of 30 dB model VAT-30 from MiniCircuits [26]. Then, this signal is low-pass filtered at 
100 MHz with the filter model SLP-100+ from MiniCircuits [27] to eliminate possible spurious signals 
out of the band of interest. Finally the baseband signal is up-converted to the frequency of interest (LI) 
using the ZX05-U432H+ mixer from MiniCircuits [28]. To do this, the mixer needs a local oscillator 
at 1,575.42 MHz — 80 MHz = 1,495.42 MHz, which has been implemented using a surface mount 
frequency synthesizer model FSW80150-10 from Synergy Microwave Corporation [29]. This device 
is programmed at the beginning through the Spartan-3 FPGA. An analog reference clock of 10 MHz 
and 1 Vpp generated in the module Master clock is necessary in the synthesizer to feed the internal 
PLL. Figure 15 shows an acquisition with the spectrum analyzer where it is possible appreciate the PRN 
bandwidth in function of the SR parameter. The higher the SR, the higher the PRN bandwidth. 

Figure 15. Spectrum analyzer acquisitions of different PRNs generated with different SRs. 
Central frequency of 1,575.42 MHz, 2 MHz/div, span = 20 MHz and 10 dB/div. (a) SR = 1, 
and (b) SR = 5 (from [6]). 
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3.2.6.2.1. Selection Circuitry 

Once the two different noise sources (thermal and PRN) have been independently generated, these 
are injected into a selector circuitry (Figure 13). The function of this part is to select between the 
Noise Source or the PRN generator, and to provide two power levels so as to perform differential 
measurements. Figure 16 shows the hardware implementation. The first element (switch 1) is an 
absorptive mechanic switch model MSP2TA-18XL from MiniCircuits [30] to select between the two 
sources. It has an isolation of 100 dB at the frequency of operation. The selected source is divided into 
two branches through a power combiner model ZAPD-2-21-3W from MiniCircuits [31]. One of the 
branches is attenuated 0 dB and the other one is attenuated 3 dB. Finally a second absorptive mechanic 
switch (switch 2) selects the required attenuated output. 



Figure 16. Hardware implementation of the selection circuitry from [6]. 




3.2.6.2.2. Correlated Noise Unit Integration 

Figure 17 shows the metallic box where the correlated noise sources and the selection circuitry have 
been integrated. 



Figure 17. Hardware implementation of the correlated noise unit from [6]. 
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In order to minimize the area these modules have been integrated in different layers. In the bottom 
layer the Noise Source module and the selection circuitry have been implemented and the PRN sequences 
module in the top layer. The FPGA 2 synthesizes the PRN sequences. It is controlled with a DB15 
connector through the internal PC as shown in Figure 17. Three commands are required to fully control 
the correlated noise unit (Figure 13): one to select the attenuation (0 or 3 dB), another one to select the 
correlated noise (Noise Source or PRN), and the last one to select the SR of the PRN. 

3.2.7. General Description of the Mobile Unit 

Before starting with the design of the mobile unit it has been necessary to define a set of 
specifications in order to establish the necessary requirements to transport the PAU-SA and the 
Multi-frequency Experimental Radiometer With Interference Tracking For Experiments Over Land And 
Littoral (MERITXELL) [32] instruments. In order to transport the instruments, it has been necessary to 
choose the mobile unit. Since the two instruments have a considerable mass and volume, a NISSAN 
ATLEON 8.19.3 truck with a maximum weight of eight tons was selected. The elevator tower for 
measurement purposes is eight meters height and it has azimuth and elevation movements 0° ^ 9 < 150° 
and —180° ^ (ft ^ +135°. It is compatible to work with both instruments, but only one is able to 
be operated while the other one is parked. The elevator tower has four positions: up or measuring, 
down or parked, calibration or looking to the internal absorber wall, and change the radiometer. All 
these commands are sent through the external computer, and finally controlled via a Programmable 
Logic Controller (PLC) located in the control panel in the truck. Moreover, the mobile unit has 
four stabilization legs manually controlled covering the maximum surface allowing to work with an 
instrument at eight meters high, under wind speed conditions up to 100 km/h. Both the elevator tower and 
the stabilization legs work with a hydraulic unit. In addition to providing the truck with an elevator tower 
for measurement purposes, it has an enclosure to storage and transport the instruments. A microwave 
absorber wall has been placed inside the mobile unit for calibration purposes, one for each instrument. 

4. PAU-SA 's Processing Implementation 

The Digital Correlation Unit (DCU) computes all possible cross-correlations between antenna pairs at 
each polarization by counting the number of samples with the same sign (1 bit), obtaining the so-called 
correlation counts matrices (N Cm n ). Figure 11 shows a single polarization correlation counts matrix 
structure: the upper triangle part of the matrix contains the I/I correlations between pairs of signals (real 
part of the visibilities samples), while the lower triangular part contains the I/Q correlations (imaginary 
part of the visibilities samples). The diagonal contains the cross-correlations of the I and Q components 
of the signals from the same antenna. These cross-correlated are related in MIRAS to the quadrature 
errors, but in PAU-SA, since the I/Q demodulation is performed digitally, there are no quadrature errors, 
and the diagonal is always zero. This 25 x 25 matrix is then completed by adding a right column 
containing the cross-correlation between the in-phase component and 0's (I-'0'), and the bottom row with 
the cross-correlation between the quadrature component and 0's (Q-'0') (both for calibration purposes). 
These values are used to compensate for threshold errors in the comparators. The total number of samples 
Ncmax is added at the right bottom element, which is used to normalize the whole N Cm n matrix. It can 



Sensors 2012, 12 



7759 



be observed that the number of correlation counts is an integer, with respect the equation: 0 < N Cmtn < 
N Cmax . Recalling that the sampling frequency (F s ) is 5.745 MHz, and that the maximum integration 
time (Ti nt ) is 1 s, then the maximum number of counts or samples (N Cmax ) is: 

N Cmax = F s • T int = 5.745 M Samples (26) 

Finally, these matrices, and the power measurements are sent to an external PC in real time, where 
the calibration and the image reconstruction algorithms are implemented to retrieve the brightness 
temperature image. 

The first step consists of computing the cross-correlation (c) by normalizing the correlation count 
matrix N Cm ^ n to the maximum number of possible samples (N Cmax ), at each polarization (V, H and V/H): 

N c 

Cmax 

where m and n indicate the corresponding receivers. Once c is computed, it is necessary to compensate 
the offset introduced in the measurements by the up-counters and apply a scaling factor to obtain 
the digital correlation (Z and \i correspond to either the real or imaginary parts of the digital and 
the normalized cross-correlation, depending if the cross-correlation is performed between the in-phase 
components, or between the in-phase and quadrature components.) (Z)\ 

Z m ,n M C m,n offset m?n ) (28) 

The digital correlation (Z) is a real number, for which it can be assumed that 0 < Z < 1. On one 
hand, when two signals are uncorrelated, the correlation c m ^ n = 0.5 and Z m ^ n = 0, and on the other 
hand, when two signals are correlated, the correlation c m ^ n = 1 and Z m ^ n = 1, therefore in the absence 
of hardware errors [33] 

Zm.n 2 I C m n — I (29) 



2 . 

The normalized cross-correlation (/x) between two Gaussian signals sampled with 1 bit (2 levels) is 
related to the digital correlation by the following expression: 

H min = sin (^Z m , n j (30) 

However, this relationship is only valid for ADCs having a zero offset (ideal case). Taking into 
account the sample threshold error, to correct the correlation offset, the normalized correlation \i is 
calculated with a non-linear relationship (Equation (31)) using Equation (30) as initial solution, and 
using an iterative method (fixed point iteration or Newton-Raphson) [34,35] until the difference between 
two iterations is smaller than, for example, 10 -6 . 

Mm,n = Sin ( — (Z mn H (fa,n^Ql + l^m.nXoi ~ 2X 0 iY 0 i)^\ ] (31) 

In Equation (31), X 0 i is calculated in the last column of the matrix of correlators (26 th column) as 
the cross-correlator between the in-phase components (/) and a signal equal to all zeros. Similarly, F 0 i 
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is computed in the last row (26 th row) as the cross-correlator between the quadrature components (Q) 
and the signal to all zeros. After the fixed-point iteration, the new normalized correlation matrices are 
reorganized to compute the 25 x 25 visibility matrices. The upper and bottom parts of the diagonal are 
combined to form the normalized complex cross-correlation. 

where the subscripts stand for the correlation between two in-phase n, and quadrature in-phase iq 
contributions. This relationship can be also be expressed as: 



V 1 SYSm 1 SYS n 



(33) 



where 3?e and 9m are the real and imaginary part operators respectively, and (0) and r^ n (0) are 
the Fringe- Wash Function (FWF) at the origin for the corresponding pair of receivers indicated by the 
sub- scripts, and the system temperatures of the denominator are given by: 

TsYS = r A + T REC (34) 

where T A is the antenna temperature including the contribution from the ohmic losses t]q (T a = t]qT a + 
(1 — Vn)T p h), and V m ^ n represents the corrected visibility: the m and the n receiver assuming that both 
receivers are at the same physical temperature T ph , which is given by: 

V r 



t+ \ (35) 

fo J 

where the overbar in the FWF means normalization at the origin, that is: 

r a @ (t) 

~ ( 4.\ m,n\ / /o/c\ 

r m ,n(t) = ^rzr (36) 

' 771,71 V U J 

Once the normalized complex cross-correlations n m ^ n have been calculated for every combination of 
antenna polarizations, they are arranged as 25 x 25 visibilities matrices. Then, the normalized visibility 
function can be derived from Equation (37). 

' V T SYSm T SYS n (37) 

Finally, the visibility samples must be corrected for phase and amplitude errors [10,36]. The main 
differences with the MIR AS instrument calibration are: 

1. Since the I/Q demodulation is performed digitally, quadrature errors are zero and do not have 
to corrected. 

2. System temperatures are measured with digital PMS, therefore they are insensitive to offset and 
slope drifts as opposed to their analog counter parts, 
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3. Visibility offsets are measured with an internal matched load and by looking to an external 
absorber, and 

4. Phase and amplitudes of the visibility samples are measured using either a centralized noise source 
or a pseudo-random noise sequence PRN [36]. 

- 1 



V m ,n — fc,n ' \/Tsy S m TsY S ,., 



De-normalization 



9m,r, 



(38) 



Phase calibration 



Amplitude calibration 

Afterwards, these visibility samples must be re-ordered and assigned to the corresponding baseline 
((u,v) point) before applying the image reconstruction algorithm. Since it has not been possible 
to measure the antenna patterns when mounted in the structure, antenna pattern mismatches are 
compensated for by using a sort of Flat Target Response as described in [10]. 

5. Inter-Comparison between MIRAS and PAU-SA 

Table 1 shows an inter-comparison table between both instruments to show up the potential 
improvements over the current MIRAS design that have been implemented and tested. 

Table 1. Comparative between MIRAS and PAU-SA. Table from [6]. 



N° Parameter 



MIRAS/SMOS 



PAU-SA 



Comments 



1 Altitude 



Global observation, Low 
Earth Orbit (LEO): orbital 
altitude of 763 km, 3 days 
equatorial revisit time 



On-ground 



2 


Frequency 


L-band (1,400-1,427 MHz) 


LI -band (1,575. 42 MHz) 


Same frequency for Radiometer and 




operation 


band is protected for passive 


GPS signal 


GNSS-Reflectrometer 






observations 






3 


Bandwidth 


19 MHz 


2.2 MHz 




4 


Number of 


4 m 


1.3 m 


Negligible spatial correlation effects 




antennas per arm 









5 Number total 
antennas 



69 



31 



8 x 3 + 1 = 25 for Radiometer, 
3 center plus 3 additional = 7 
antennas for GNSS-Reflectometer, 
3dummy antennas, 1 at the end of 
each arm 



7 Antenna type 



Patch antenna without 
dielectric substrate and V & 
H polarizations 
(non-simultaneous) 



Patch antenna without 
dielectric substrate and V 
& H polarizations 
(simultaneous) 



Full-polarimetric (non-sequential) 



8 Antenna spacing 



0.875 A at 1,400 MHz, 21 
cm wavelength 



0.816Aat 1,575.42 MHz, 
19 cm wavelength 



Increase the alias-free field of view 
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Table 1. Cont. 



N° 


Parameter 


MIRAS/SMOS 


PAU-SA 


Comments 


9 


Receiver type 


1 per element 


1 per polarization 
(2 per element) 


Full-polarimetric possible 
(non-sequential) 


1(J 


lopology 
of the LO 
down-converter 


Distributed local oscillator 
(LO) (groups of 6 elements) 


Centralized reference 
clock + Internal LO 
generator 


Elimination of correlation offsets 
due to LO noise leakage. 


1 1 


Quantization 


1 bit IF sampling depending 
upon the noise uptake level 
(Inside the LICEF ) 


8 bit IF sub-sampling 
using an external ADC 


(8 bits) for I/Q conversion and 
(1 bit) to power measurement 


12 


1/Q down- 
conversion 


Analog 


Digital 


Mass reduction, no quadrature 
errors (calibration not required) 


13 


Frequency 
response shaped 
by 


Analog RF filter 


Digital low- pass filter 


Mass reduction, quasi perfect 
matching, no temperature and 
frequency drifts 


14 


Power 

measurement 
system (PMS) 


Analog (diode detector) 


Multibit Digital (FPGA) 
Computation 


Mass reduction, no temperature 
drifts 


15 


Digital 

Correlated Unit 


Fclk = F s 


Fclk » F s 


Clock frequency (Fclk) much 
higher than sampling frequency 
(F s ) allows hardware reuse and 
compute full-polarimetric 
correlation matrices in one snapshot 
inside FPGA 


id 


Image 
capabilities 


Dual-polarization or 

full-polarimetric 

(sequential) 


Full-polarimetric 
(non- sequential) 


Necessary for GNSS-R applications 


17 


Integration time 


1.2 s 


Variable: 4 values 1 s, 
0.5 s, 100 ms, 10 ms 




18 


Correlated Noise 
Injection 


Distributed (Noise Source) 


Centralized (Noise Source, 
or PRNs) 


Using PRNs independent number of 
receivers (simpler and more flexible 
calibration) 



One of the main differences between MIRAS and PAU-SA is the frequency of operation. L-band 
radiometer should operate in the 1,400-1,427 MHz "protected" band as MIRAS. However, PAU-SA is 
an instrument concept demonstrator, and to minimize the hardware requirements the GPS reflectometer 
and the L-band radiometer share the same front-end and frequency band. Although sharing the same 
front-end provides a significant hardware reduction, there are some drawbacks. The first one concerns 
to the possible interference that GPS signals can introduce in the radiometric measurements. This 
non-optimal operation frequency for the radiometer instrument has an impact on the radiometric 
measurements, introducing some errors. This issue is analyzed in Section 5.1. The second one concerns 
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the bandwidth. In the case of MIR AS, the bandwidth is limited by the protected band with a maximum 
value of 27 MHz, and an effective noise bandwidth of ^19 MHz. For the PAU-SA instrument the 
bandwidth is 2.2 MHz imposed by the IF frequency of the GP2015 chip from Zarlink being used, and 
by the SAW filters used for the implementation of the receiver chain. This reduction in the bandwidth 
has an impact on the radiometric sensitivity that can only be compensated by increasing the integration 
time, which is not critical in a ground-based instrument. 

Each arm of the MIR AS instrument is approximately three times longer than the PAU-SA ones: ^4 m 
with 23 elements in front of 1 .3 m with just 8 elements. The total number of antennas in MIRAS is 69 and 
25 in the case of PAU-SA. This decision was taken for two reasons: the first one is due to the use of fixed 
non-foldable arms in order to simplify the mechanical complexity, and the second one was pragmatic: 
to be able to take the instrument out of the laboratory, where it was assembled. Since the bandwidth of 
PAU-SA is narrower than that used in MIRAS and the arm length is smaller, spatial decorrelation effects 
modeled by the FWF are negligible. This factor is quantified in Section 5.2. One of the improvements 
of PAU-SA is the additional dummy antenna at the end of each arm to improve the antenna pattern 
similarity (Figure 2). 

Concerning the antenna type and separation, they are quite similar in both instruments. For instance, 
both MIRAS and PAU-SA use patch antennas without dielectric substrate for the V- and H-polarizations. 
For hardware simplicity, MIRAS has sequential acquisitions since the receiver is shared for both 
polarizations. In the case of PAU-SA, each polarization has its own receiver channel. Therefore, 
continuous acquisitions at both polarizations can be obtained simultaneously, as well as mixing 
the polarizations. 

In order to increase the AF-FOV, the minimum distance between element spacing is kept to the 
minimum, only limited by the receiver's size. MIRAS has an antenna spacing of 18.75 cm, corresponding 
to 0.875 A at 1,400 MHz. In the case of PAU-SA the antenna spacing is reduced to 15.5 cm = 0.816 A at 
1,575.42 MHz, the GPS LI signal. 

Moreover, the distribution of the LO used in the down-converter is different in both instruments. 
In MIRAS the LO is fed to groups of 6 elements forming an arm section, whereas PAU-SA uses a 
centralized reference clock of 10 MHz, and the LO is generated by means of a PLL inside each receiver 
to minimize offsets coming from common LO noise leaking through the mixer. 

MIRAS 'quantification scheme uses 1 bit sampling at intermediate frequency IF [37] depending upon 
the noise uptake level inside the Light Cost Effective Front-end (LICEF). PAU-SA uses 8 bits ADC 
using IF sub-sampling techniques to down-convert and demodulate simultaneously. Using 8 bits, it is 
possible to filter, demodulate, and perform the power estimation, and finally use only 1 bit to obtain the 
three complex correlation matrices V, H, V/H. 

Due to the large number of receivers in interferometric radiometers, it is advisable to obtain a quasi- 
perfect matching of the frequency responses, mass reduction, and to eliminate temperature and frequency 
drifts as much as possible. For these reasons, the most important contributions in PAU-SA are focused 
on the replacement of analog by digital subsystems being the most important: 

• I/Q down-conversion to eliminate quadrature errors, 

• Digital filtering, replacing the narrow RF filter by a digital IF filter, to obtain a mass reduction, a 
quasi perfect matching, and eliminating thermal and frequency drifts, and 
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• Digital Power estimation, eliminating the classical Schottky or tunnel diodes to achieve a mass 
reduction, and eliminating temperature drifts and aging. 

All these subsystems and the DCU that computes the full cross-correlation matrix (V, H and V/H) 
have been implemented in a Virtex 4 FPGA. In this case, since the clock frequency is much higher than 
the sampling frequency, hardware reuse techniques are used to compute the full-polarization matrices in 
each snapshot. 

The imaging capabilities can be sequential dual-polarization or full-polarimetric in MIRAS, and 
non- sequential full-polarimetric in PAU-SA. The use of both polarizations simultaneously is also 
necessary to compose the reflected LHCP GPS signal. 

The integration time is fixed for MIRAS with a value of 1.2 s, while PAU-SA has predefined 4 values: 
10 ms, 100 ms, 0.5 s and 1 s for test purposes. 

Finally, MIRAS uses a classical correlated noise injection method for calibration proposes. Due 
to the large number of receivers to feed, it is necessary to use several noise sources distributed along 
the instrument increasing the hardware complexity, and introducing additional noise. In addition, the 
distributed noise injection is not capable of calibrating all sorts of errors. To overcome this problem 
PRN signals can be used instead of a centralized noise source for calibration purposes. Moreover, 
since the PRN signals are deterministic and known, new calibration approaches are feasible through 
the correlation of the output signals with a local replica of the PRN signal, leading to the estimation 
of the receivers' frequency responses and the FWF [36]. PAU-SA has the possibility of using both the 
classical noise-injection method or this new technique with PRNs. 

5.7. Impact of the Frequency Operation on the Radiometer Part 

The receiver's operating frequency is defined by the LI signal of the GPS signal (1,575.42 MHz), 
which is also suitable for SSS estimation. On one hand the GNSS signals use spread spectrum techniques. 
After the scattering on the sea surface, the power of these signal is at least 23 dB below the thermal 
noise. For this reason, and thanks to the 30.1 dB correlation gain, GNSS-Reflectrometry can detect the 
GPS signal when the correct C/A code is applied. On the other hand, from the radiometer point of 
view, the noise signal to be detected is at least 23 dB above the GPS signal so the radiometric error 
induced is minimum, and it only occurs in the directions of specular reflection which are known a 
prior. Therefore, it is possible that both the radiometer and the reflectometer share the same receiver. 
Although the spread- spectrum of the GPS modulation is at least 23 dB below the noise level (PAU-SA 
as a radiometer), its impact has to be quantified specially since the synthesized beam is very narrow. The 
GPS signal includes the P, C/A and M codes. Figure 18(a) shows the auto-correlation function (ACF) of 
P and C/A codes. 

The first one is the P code, having the largest repetition period of 6.1871 x 10 12 bits long 
(6,187,100,000,000 bits, ^773.39 gigabytes) repeating once a week, and achieving a ^75 dB 
compression gain with a 22 MHz band, and the received power level is ~— 133 dBm (—143 dBm in the 
~2.2 MHz C/A code bandwidth) [39]. The C/A code has a period of 1,023 bits, repeating every 1 ms, and 
achieving a 30 dB of compression gain with a 2.2 MHz, and the received power level is ~— 130 dBm. M 
code is still experimental and, in addition, it distributes its power at the edges of the band (Figure 18(a)), 
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having less effect even than the P code. From the point of view of the radiometric measurements, the 
correlation gain does not have any effect, and only matters the power density associated at each code. 
For this reason the GPS signal to be taken into account is the C/A code, being 10 dB higher than the 
P code. To measure the impact of the GPS signal in the PAU-SA's radiometric measurements, the 
PAU-SA's antenna pattern presented in the previous chapter has been analyzed in two extreme cases: 
the worst case is when the specular reflection point is at the antenna's boresight (maximum directivity 
of the synthetic antenna pattern: 46.5 dB), and the best situation comes from the interference at the 
edge of the AF-FOV (minimum directivity of the synthetic antenna pattern: 45.5 dB), only 1 dB below 
the maximum. Figure 18(b) shows the contribution of the GPS signal (C/A code) to the radiometric 
measurements AT considering a sea surface reflection coefficient of T = 0.7, and the PAU-SA system 
parameters. As it can be appreciated, independently of where the interference comes from, the GPS 
signal has a high contribution to the radiometric measurements (between 170 K and 220 K respectively). 
This means that, in the case of PAU-SA, to perform radiometric measurements the instrument must 
be pointed to the North, where the density of GPS satellites is minimum due to their orbital plane 
distribution. Moreover, if a GPS satellite interferes in the radiometric measurements, it will appear as a 
"point source" in the retrieved image, which could be subtracted by making measurements at different 
time. (GPS satellites will appear at a different position). 

Figure 18. (a) LI ACF of the P code (in blue), and the C/A code (in red) from [38], 
(b) estimated AT B contribution in the systematic T B image due to GPS C/A code in the 
PAU-SA's radiometric part as a function of the GPS position in the FOV from [6]. 
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5.2. Impact of the Spatial Decorrelation Effects in the Visibility Function 



The Fringe- Wash Function FWF accounts for the spatial decorrelation of the signals coming from 
a given direction at a given baseline. Precisely, the FWF is related to the frequency responses of the 
receivers forming the baseline. 
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The amplitude of the FWF can be modeled around the origin by means of a sine function as: 

\G kj {r)\ « A • sinc(B • (r - C)) (39) 

where ^4 • sinc(B • (7) is the amplitude at r = 0, B is the noise bandwidth, and (7 is the value of r in 
which the fringe- washing function is maximum. The transit time r from given (£, 77) directions to the 
antennas forming the baseline (u, v) is 

T = -<±^1 (40) 



/o 

The maximum value is found for the largest baseline Au max , and the (£, 77) direction distance farthest 
away from the boresight within the AF-FOV, | (£, 77) | max = 0.519 (Figure 19(a)). 

Figure 19. (a) PAU-SA's AF-FOV showing the maximum distance in the plane (£, 77), and 
(b) FWF comparison between MIRAS and PAU-SA instruments, (from [6]). 
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The maximum antenna separation is given by: 

Au max = 2V3N EL d (41) 

where N EL is the number of antennas in each arm and d is the distance between adjacent antennas. In 
the case of PAU-SA N EL = 8, d = 0.816 A, and Au max = 22.61 A. Figure 19(b) shows a comparison 
between the FWF of MIRAS and PAU-SA systems. In this case, since the bandwidth in MIRAS is 
about 10 times larger than PAU-SA, and the array size is approximately 3 times larger than PAU-SA, the 
FWF is ^30 times narrower in MIRAS than in PAU-SA. As it can be appreciated in Figure 19(b), the 
PAU-SA's FWF in the evaluated range is very close to one (FWF ~ 0.9995), so that this parameter is 
negligible in the PAU-SA system. 

6. Instrument Characterization and Experimental Results 

This section is divided into three parts. The first part presents the temperature characterization of the 
instrument, the second part is the characterization at baseline level in an anechoic chamber, and the last 
part is an outdoor experiment with the instrument. 
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6.1. Thermal Control Characterization 

Since receivers exhibit phase and amplitude drifts due to temperature changes, it is necessary to 
stabilize and control the temperature of the instrument. The better the temperature control, the longer 
the inter-calibration period will be. The PAU-SA's stability temperature is controlled by Proportional 
Integral Derivative (PID) devices located in each arm, and in the central part or hub. Moreover, it has an 
air control by means of fans to distribute the air along the instrument and forcing the air circulation as 
much as possible, but in the case of the HUB it is not sufficient, being the warmest the central elements. 
Figure 20 shows a plot of the control temperature during one day long in a measurement campaign. As 
it can be appreciated, in the stable transition there is a temperature excursion of about 10 °C, however, 
the maximum standard deviation of each sensor which is acceptable is only 0.27 °C. 

Figure 20. Sample time evolution of the PAU-SA's receivers acquired during the day 
(2011-02-15) from [6]. 
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6.2. Measurements at Baseline Level 

The tests performed at baseline level are: the characterization of the radiometer noise and its stability 
through the 'Allan's variance" [40], and the characterization of the radiometer resolution. 

6.2.1. Radiometer Stability 

In order to determine the optimum range of integration times for best use of the system, the 
characterization of the radiometer noise and stability has been performed measuring the Allan's 
variance [40] given by: 

° 2 A(r)= l -((ji n+1 -Ji n ) 2 ) (42) 

where r is the integration time, and JI n is the n th fractional frequency average over the observation period. 
This method consists of the determination of the Allan's variance versus integration time allowing to 
determine the different types of fluctuations of the radiometer output signal. In particular, the range of 
integration times for the optimum setup of the calibration and measurement is shown in Figure 21. 
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Figure 21. (a) Normalized power variance versus number of samples, (b) normalized 
correlation (real and imaginary parts) vs. number of samples. Pictures from [6]. 
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It shows the evolution of the variance of both channels versus the number of samples. The ideal case 
(thermal Gaussian noise) must decrease monotonically decreasing function as: 

<* = idb <43) 

However, due to system instabilities, after a given value of N, the variance o\ increases again. 
Taking into account that a sample represents a measurement of 1 s, the optimum number of samples 
is the value of minimum variance and determines the maximum integration time that may be used, 
without degradation due to radiometer drifts. On one hand, Figure 21(a) plots the power variance vs. 
the number of samples. The maximum integration time for calibration purposes using correlated noise, 
and measuring the visibilities de-normalization procedures is 15 s. One the other hand, Figure 21(b) 
shows the normalized correlation vs. the number of samples to determine for phase and the ADC offset 
calibration. This time is ~200 s. Therefore, gain fluctuations are dominant and are the ones that limit 
the maximum integration time. 

6.3. Radiometer Resolution Validation 

The characterization of the radiometric resolution has been performed with the set up shown in 
Figure 22(a). A noise source or a matched load are connected to the input of a non-resistive power 
splitter producing correlated noise. Two adjustable attenuators are connected to the power splitter outputs 
producing uncorrected noise and attenuating the correlated noise generated by the matched load. The 
complex correlation is then measured for different values of the attenuation and for different phases. 
Phases have been randomly obtained by turning OFF and ON rapidly the PLL, so each time the receiver's 
PLL locks to a different phase. 

Concerning the shape of the circle is clearly distinguished up to 2 x 20 dB attenuators, and since the 
I/Q demodulation is performed digitally, there are neither quadrature errors, nor amplitude unbalances 
between branches that need to be corrected [41]. The radius of the circles can be determined estimating 
the correlation /i of two digital signals (Equation (44)): 

» = ^ = v 3r (44) 

J- SYS 1 A^~ 1 REC 



Sensors 2012, 12 



7769 



where /i is equal to the corresponding analog signal power (T c ) normalized to the system temperature 
T S YS- Neglecting the outgoing noise waves from each receiver front-end that couple into the other 
receiver, T S ys can be estimated as (Equation (45)): 

T r „ / 1 



SYS 



2L 



T ph [ 1 



L 



REC 



(45) 



where L is the attenuation introduced, T ph is the physical temperature, and T c = 290 K (matched 
load). Evaluating the previous equations with L = 20 dB, T ph = 290 K, and T REC ~ 250 K than, 
Tsys == 538.55 K, and the correspondent /i = 2.7 • 10" 3 , which is very close to the measured value 
(Figure 22(b)). 

Figure 22. (a) Sensitivity measurement test set-up, and (b) measurement sensitivity circle 
with a matched load and an attenuation of 20 dB with an integration time of 1 s. Pictures 
from [6]. 
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Moreover other tests have been performed in the anechoic chamber [42] to measure the normalized 
visibility by scanning the pairs of receivers from —90° to 90°, at the polarization from V to H as shown 
in Figure 23(a). Results are contrasted with theoretical results through Equation (4). Considering that 
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the fringe- wash function is negligible (1/B <C maximum transit time) [43] and that both antenna patterns 
can be considered to be equal, the normalized visibility can be determined by: 



V(d) 



4 7rr 2 4^ 



D R t(6)e 



-j27rf sin^ 



(46) 



In this situation, the theoretical real and imaginary normalized correlations should vary according to: 
real part in Equation (47) and the imaginary part in Equation (48) as shown in Figure 23(b). 



Hi 



2nd . 
sin | — ^— sin v Q 

[2nd 
cos(— 



sin 6c 



(47) 
(48) 



Due to the previous considerations, in Equation (46) the absolute value of the normalized visibility is 
the antenna radiation voltage pattern t(9) when the transmitting antenna has an orthogonal polarization. 
Figure 23(c) shows that the normalized correlation sweeps the emitter with different polarization, in the 
center with opposite polarization. 

Figure 23. (a) Picture of the baseline level measurements during anechoic chamber test, 
and normalized correlation (real and imaginary) measurements with baseline rotating at 
(b) X-axis, and (c) Z-axis. Pictures from [6]. 
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Figure 23. Cont. 
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6 A. Experimental Results 



This section summarizes a couple of experimental results carried out with the PAU-SA instrument, 
calibration procedures and extended results are presented in detail in [10]. The instrument performance 
has been evaluated in terms of AF-FOV, angular resolution, radiometric resolution, radiometric 
precision, and radiometric bias. Moreover, since PAU-SA operates at the GPS LI band, maps of the 
GPS satellites constellation have been obtained. To check at first glance the instrument's health, firstly 
a PRN source has been transmitted in order to confirm that the source was imaged in the right direction 
and also determine the AF-FOV. 

Figure 24. Measurement setup using a PRN signal as point sources at 10 m of the instrument, 
(a) single transmission moving the instrument in azimuth and elevations angles to determine 
the AF-FOV, and (b) recovered image moving the instrument 20° in azimuth an applying 
near-field to far-field compensation (from [6]). 
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Figure 24(a) shows the PAU-SA instrument in the top of the 8 m height robotic arm and the beacon on 
the ground. Since it is easier to move the instrument than the beacon, the PAU-SA instrument was moved 
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±10° and ±20° both in azimuth and elevation. Figure 24(b) shows one of the images retrieved when the 
instrument was moved 20° in azimuth, which can be compared with the theoretical one (Equation (49)), 
and determines the AF-FOV to be approximately around 48°. 

- • AF - FOV = arcsin ( -jt- - 1 J = 24° (49) 
2 \\/3d J 

The second test has been the determination of the angular resolution transmitting two point sources 
to the instrument and using a hexagonal inverse Fourier transform, and a rectangular window with 
an angular resolution of A£ = 0.1 corresponding with 5.7°. The radiometric specifications have been 
achieved considering the same antenna pattern and neglecting the FWF and the Corbella's term: 

• The radiometric resolution is the standard deviation of the time fluctuation of a given observable. 
It is the minimum change detectable by the instrument and it is computed as in Equation (50). It 
has been found to be a Vjh =1.9 K at both polarizations. 

• The radiometric precision is the systematic error in each pixel. An average radiometric precision 
value is computed for the whole alias-free field-of-view image as the RMS value of the brightness 
temperature computed from the average of 80 snapshots of 3 s integration time each (total 240 s), 
so as to achieve a negligible radiometric resolution). It is estimated by using Equation (51) a v = 1 .2 
K and a h = 2.0 K, and 

• The radiometric bias is the spatial mean of the computed brightness temperature minus a reference 
temperature determined by Equation (52). It is found to be —1.6 K and —1.8 K, at vertical and 
horizontal polarizations, respectively. 



\ T ~ M / ££i (Tb&kU) - (f B (£,ri,t)) t ) 2 

resolution ~ \/ ^ W U J 



, T {(TB^i,rh,t)) t -T B (^rH)) 2 

t-^-L precision ~ y jV 1 

l M ~ 

AT bias = jjJ2 « Tb ^ *>, U))t ~ T B (£, V )) (52) 

2=1 

Finally, since the instrument was conceived as a technology demonstrator, commercial GPS chips 
operating at LI band were used, making it possible to track the GPS constellation. Figure 25(a) shows 
the instrument pointing to the sky seeking satellites, Figure 25(b) shows the map of the expected GPS 
satellites' paths in this time and location, and Figure 25(c) shows a video animation of the satellites' 
movements appearing as points sources and following the expected tracks. 
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Figure 25. (a) PAU-SA pointing to the GPS satellites, (b) Map of GPS satellites' paths 
as seen from the test location on March 30th, 2011. Sequential of images recovered by 
PAU-SA every 44 min: (c) animation (click on the image) of the GPS satellites' movements 
UTC 12:00:03. Note: £ = 0, r] > 0 corresponds to the geographic north from [6]). 
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7. Conclusions 

This paper has described the PAU-SA instrument. It is a synthetic aperture radiometer that has been 
designed and tested to study potential improvements in SMOS follow-on or in other missions using 
synthetic aperture radiometers as GAS [44], PATH [45] or GIMS [46]. A comparison table between 
MIRAS and PAU-SA has been presented, describing the main novelties of the PAU-SA instrument 
design. These contributions have been focused on the replacement of analog by digital subsystems 
such as: I/Q down-conversion, digital filtering, full-correlation matrix (V, H and VH), power estimation 
implemented in an FPGA, and the use of both a centralized noise source, and PRN signals for calibration 
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purposes. This last one has been one of the most remarkable contributions in the hardware design, 
since with this method it is possible to feed a large number of receivers using a centralized topology. 
Once the instrument has been assembled and calibrated, imaging results can be obtained. Instrument 
characterization has focused on the measurement of Allan's variance, the angular resolution, and the 
radiometric performance. Finally, imaging results recovery of point sources using PRN sequences and 
imaging real GPS satellites have been performed. However, the radiometric performances achieved so 
far are not good enough for the sea surface salinity retrievals originally foreseen. 
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